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Abstract. Two methods to change a quantum harmonic oscillator frequency without 
transitions in a finite time are described and compared. The first method, a 
transitionless-tracking algorithm, makes use of a generalized harmonic oscillator and a 
non-local potential. The second method, based on engineering an invariant of motion, 
only modifies the harmonic frequency in time, keeping the potential local at all times. 
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1. Introduction 

Changing the external parameters of the Hamiltonian is a fundamental and standard 
operation to probe, control, or prepare a quantum system. In many cases it is desirable 
to go from an initial parameter configuration to a final one without inducing transitions, 
as in the expansions performed in fountain clocks pQ. In fact most of the current 
experiments with cold atoms are based on a cooling stage and then an adiabatic drive 
of the system to some desired final trap or regime [2]. These "transitionless" [3J, or 
"frictionless" [I] adiabatic processes may require exceedingly large times and become 
impractical, even impossible [2], or quite simply a faster process is desirable, e.g. to 
increase the repetition rate of a cycle, or a signal-to-noise ratio. This motivates the 
generic objective of achieving the same final state as the slow adiabatic processes, 
possibly up to phase factors, but in a much shorter time. One may try to fulfill that 
goal in two different ways: (a) designing appropriate "parameter trajectories" of the 
Hamiltonian from the initial to the final times, or (b) applying entirely new interactions 
that modify the Hamiltonian beyond a simple parameter evolution of the original form, 
for example by adding different terms to it. In this paper we shall analyze and discuss, 
for the harmonic oscillator, two recently proposed methods whose relation had not been 
investigated. It turns out that they actually implement these two different routes. While 
most of the treatment is applicable to an "abstract" harmonic oscillator, we shall discuss 
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physical implementations specific of ultracold atoms or ions. Indeed, harmonic traps 
and their manipulation are basic working horses of this field. 

For the harmonic oscillator the parameter we consider is the trap frequency, 
which should go from w to in a time tf, preserving the populations of the levels, 
P n (tf) = -Pn(O). "n" labels the instantaneous n-th eigenstate of the initial and final 
harmonic oscillator Hamiltonians, 

#o(0)|n(0)) =^ (n + l/2)|n(0)), 

H {t f )\n(t f )) = hu f (n + l/2)\n(tf)). (1) 

One of the methods we shall discuss here relies on a general framework set by Kato 
in a proof of the adiabatic theorem [5], and has been formulated recently by Berry 
[3]. We shall term it "transitionless-tracking" approach, or TT for short; the other one 
P, [7] engineers the Lewis-Riesenfeld invariant [8] by an inverse method [9] to satisfy 
the desired boundary conditions; we shall call this method "inverse-invariant", or II 
for short. In the basic version of TT the dynamics is set to follow at all intermediate 
times the adiabatic path defined by an auxiliary Hamiltonian H (t) (in our case a 
regular harmonic oscillator with frequency u(t) and boundary conditions uj(0) = ujq and 
ujf = u(t = tf)), and its instantaneous eigenvectors \n(t)), up to phase factors. Instead, 
in the II approach the auxiliary object is an engineered Lewis-Riesenfeld invariant I(t) 
set to commute with Hq(0) at t — and with H^itf) at tf. In both cases intermediate 
states may be highly non-adiabatic with respect to the instantaneous eigenstates of the 
Hamiltonians actually applied, Htt^) and Hu(t). 

We shall provide first the equations characterizing the two approaches and then 
comment on possible physical implementations. 



2. Transitionless tracking algorithm 

2.1. General formalism 

For the general formalism we follow |3J closely. Assume a time-dependent Hamiltonian 
H (t) with initial and final values (CD), instantaneous eigenvectors \n(t)) and eigenvalues 
E n (t), 

H (t)\n(t)) = E n (t)\n(t)). (2) 

A slow change would preserve the eigenvalue and eigenvector along the dynamical 
evolution times a phase factor, 

|Vn(*)> = exp j-1 J*dt'E n (t') - ^dtf(n{t')\d t ,n{t'))} W)). (3) 

We now seek a Hamiltonian H(t) such that the adiabatic approximation \ip n (t)) 
represents the exact dynamics, 

m\Mt)) = H{t)\i/> n (t)). (4) 
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H(t) (which is Htt if distinction with the other method is needed) is related to the 
corresponding unitary operator by 

iWtUit) = H(t)U(t), (5) 

H(t) =i%{d t U{t))U\t). (6) 

Choosing 

U(t) = X>p {-^ f Q dt'E n {t') - jf* dt> (n(t')\d t/ n(t'))} Ht))(n(0)\, (7) 
we find from (J6]), 

H(t) = J2 \n)E n {n\ + i%J2(\d t n) (n\ - (n\d t n)\n)(n\) = H + Hi, (8) 

n n 

where we have simplified the notation, \n) = \n(t)). It is also possible to choose other 
phases in ([3]) [3J. The simplest case is U(t) = l n (^))( n (0)|, without phase factors, 
which leads to H(t) = ihJ2 \d t n)(n\. Note that with this choice H (t) has been formally 
suppressed in H(t) but still plays a role through its eigenf unctions \n(t)). 

2.2. Application to the harmonic oscillator 

We now apply the above to the harmonic oscillator 

H (t) = fjlm + uj(tfx 2 /2m = hu{t)(d\d t + 1/2), (9) 

where a t and df are the (Schrodinger picture!) annihilation and creation operators at 
time t, 



£ = j 2mW) i al + at) ' (10) 
P =i\h^^ L {a\-a t ) ) (11) 



mu(t) ( i \ , . 



This time dependence may be misleading and a bit unusual at first so we insist: since 
the frequency depends on time the "instantaneous" ladder operators d t , d\ create or 
annihilate different "instantaneous" states, adapted to the corresponding frequency 
Thus, ladder operators with different time labels do not commute in general, although 
some combinations, e.g. those equivalent to powers of x and/or p, do commute, as we 
shall see later. 

The instantaneous eigenstates \n(t)) can be written in coordinate representation as 
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and their derivative with respect to t is 
(x\d t n(t)) 



1 



muj(t) 2 



~ '»> + \h^^Mn - 1). (15) 



2h ~ J u(ty V 2h u(t) 

We find, using the recursion relation of Hermite polynomials and their orthogonality. 



(k\d t n) 



(n + l)(n + 2) 



w(t) 







k = n - 2 



fc = n + 2 



(otherwise) 



(16) 



so that H\{t) can be written as 



H\{t) =ihJ2 \dtn)(n\ = ^-tttXI 



1 



muj(t) „ 2 



2ft 



x |w)(n| 



+ 



' muj{t) 

2h 



xy/n\n — 1) (n| 



;i7) 



Using a< = 5Z n \/n|n — l(t))(n(t)|, and the relations between x, p, at and a\ written 
above, 



ih Cj 



l 

4 ~ 

1 UJ 



2h + 



'moj(t) 



2h 



-xa t 



= -; tt rr%P- (18) 

4 uj(t) 2uj(t) y K ' 

Using [x,p] = ih, we finally write the Hamiltonian Hi(t) in the following simple forms 
H x (t) = ~(xp + px) = th^[a 2 - (a 1 ") 2 ]. (19) 



In the last expression the subscript t in a and a^ has been dropped because the squeezing 



combination a 2 



^2 



is actually independent of time, so one may evaluate it at any 



convenient time, e.g. at t — 0. The connection with squeezing operators is worked out 
in the appendix. 

Hi is therefore a non-local operator, and does not have the form of a regular 
harmonic oscillator potential with an x 2 term. Nevertheless the final Hamiltonian 
H = Hq + Hi is still quadratic in x and p, so it may be considered a generalised 
harmonic oscillator |10|. 



2.3. Physical realization 

The nonlocality of Hi, with a constant prefactor, can be realized in a laboratory by 
means of 2-photon Raman transitions for trapped ions [HJ[T2]. Since we have to evaluate 
as well the possibility of making the prefactor in Hi time dependent we need to provide 
the derivation with some detail, first for a time- independent uj. 
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2.3.1. Raman two-photon transition in a trapped ion Let us consider a harmonically 
trapped two-level system in ID driven by two different lasers (with coupling strengths 
Qj and frequencies uij, j = 1,2), see Fig. [Hand Refs. [UJEE3]. The time dependent 
"Raman" Hamiltonian in the Schrodinger picture will be given by 

H R (t) = H T + H A + H inU (20) 

with "trap" (T), "atomic" (A), and interaction (int) terms 

H T = hioa^a, (21) 

H A =?iu e \e)(e\, (22) 



2 



Hm = E h ty cos ~~ k i x + fa) (\9)( e \ + \ e )(9\), (23) 

3=1 

where hu e is the energy of the excited state |e) and k, = kjit the wavevector of each laser 
which are assumed to be pointing along the principal trap direction, the x-direction. 

2.3.2. Interaction picture Let us now write the above Hamiltonian in an interaction 
picture defined by the Hamiltonian h = Ht + ?uvL\e}(e\, where l>l = (ui + U2)/2 has 
been introduced. The interaction Hamiltonian Hj = e lhot l % [H R — ho)e~ %hot l h reads 

Hi(t) = - KA\e)(e\ 

+ j2 ™k ^[«W+ at (*)] e ^K-^)* e -^| e )^| + H.c) , (24) 

j=l 2 V / 

where A = ujl — cu e , and now a(t) = ae~ lujt , a){t) = a^e tuJt are the time dependent 
Heisenberg annihilation and creation operators respectively. Note also that fast 
oscillating off- resonant e ± l ("j+^)* terms have been neglected in the rotating wave 
approximation (RWA). The parameter rjj = kjXo is known as the Lamb-Dicke (LD) 



parameter, where Xq = \jhj2muj is the extension (square root of the variance) of the 
ion's ground state, i. e., x = x (a + a)). 

2. 3. 3. Adiabatic elimination and effective Hamiltonian For a general wavefunction (in 
the corresponding interaction picture) such as 

oo 

\Mt)) = E [9n{t)\g,n) + e n (t)\e,n)} (25) 

n=0 

the differential equations of motion for the probability amplitudes g n (t) and e n (t) are 
obtained from the Schrodinger equation ihd t \ipi(t)) = Hi\ipi(t)), 

1 2 oo 

Wn(t) = ^EE n^ t+ ^\n\e-^[^V} \n')e n ,(t), (26) 

2 j=l n '=0 

i 2 oo 

ie n {t) = - Ae n (t) + - E E n J e^ 9 M\n\eM*®+*W] \ n ')g n ,{t), (27) 

^ j=l n '=0 

where Q 3 - = Uj — ljl. For large detunings, i. e., for |A| ^> Qj,u, see Fig. [H and for an 
ion initially in the ground state one may assume that the excited state |e) is scarcely 
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Figure 1. Schematic electronic and vibrational level structure for a two-photon 
transition in an ion trapped with frequency lo. uj\ and ll>2 are the laser frequencies, 
and w e the transition frequency between ground and excited states. See the text for 
further details. 



populated and it may be adiabatically eliminated. Then, setting e(t) = 0, e n (t) may be 
written as a function of the g n >{t) from Eq. ( 1271) . and substituting this result into (1261) 
there results a differential equation for the ground state probability amplitude, 

Wn(t) = Sg n (t) + — ^n,n'(t)9n'(t), (28) 
1 n'=0 

where 

s = 1 - 2 , 29 

^n,n'(i) = (n|e-^[ a( ' )+at(t) ] |n')e iM) + (n|e^ a(t)+atW ] |n')e- iM , (30) 
and where the effective two-photon Raman parameters, denoted by tildes, are given by 



5 


= UJi — uo 2 


V 


= Vl~ V2, 




= 01-02 


h 




~2 


4A ' 



(31) 

The equation for the ground state probability amplitude corresponds to an effective 
Hamiltonian 

H eff = hs\g)(g\ + ^ ^m) + aHt)] e -i(k+4>) + R ^ \ g) ^ (32) 

Note that the Stark-Shift produced by off resonant driving is included in s, which is a 
constant of motion and produces no effect on the Raman coupling between sidebands. 
We have thus adiabatically eliminated the excited state |e) ending with a Hamiltonian 
of the same form as (1241) where the transitions between electronic levels are not present. 
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2.3.4- Two-photon Jaynes-Cummings Hamiltonian in the Raman Scheme: Vibrational 
RWA Using the Baker- Campbell-Hausdorff (BCH) identity, the exponential in the 
effective Hamiltonian ( 1321) may be expanded in power series of fj (T5j 02] , 



ft - hQ 

Heff ~ — 



:~\n+n' 

n\n'\ 



(33) 



If the effective detuning is 5 = u)\ — uj 2 = 2u, the second blue sideband becomes 
resonant, and we may neglect rapidly oscillating terms in a second or vibrational RWA 
[16J. The above Hamiltonian is then simplified to a two-photon Jaynes-Cummings-like 
Hamiltonian without electronic transitions. To leading order in fj it takes the form 

where for the last step a relative phase between the applied fields = 0i — 2 = — vr/2 
has been assumed. 

2.3.5. Validity for time- dependent u Unfortunately the above formal manipulations 
and approximations cannot be carried out in general for a time dependent u. The 
interaction picture performed in I2.3.2[ in particular, assumes a constant ho- A time 
dependent one would require a more complex approach with time-ordering operators 
[T7] . Similarly, the vibrational rotating wave approximation requires the stability of the 
frequency for times larger than a period to avoid off-resonant couplings. One may still 
obtain ([34]) for a sufficiently slowly varying u, the criterion being that the change of the 
time-dependent trapping frequency in one time period T has to be much smaller than 
the frequency itself. We can write this condition as u>(t)T u(t) or 

u(t) 2 

which turns out to be the adiabaticity condition for the harmonic oscillator. Of course, 
if satisfied, the whole enterprise of applying the TT method would be useless. These 
arguments are far from constituting a proof that the TT method cannot be implemented 
for the harmonic oscillator. They simply leave this as an open question. 

3. Engineering the Lewis-Riesenfeld invariant 

In this section we describe a different method for transitionless dynamics of the harmonic 
oscillator [6j. A harmonic oscillator such as H (t) in Eq. (Q has the following time 
dependent invariant [8] 

Ht) = i gm^ + i* 2 ) , (36) 

where n = b{t)p — mbx plays the role of a momentum conjugate to x/b, the dots are 
derivatives with respect to time, and Uq is in principle an arbitrary constant. The 
scaling, dimensionless function b = b(t) satisfies the subsidiary condition 

b + u{t) 2 b = u^/b 3 , (37) 
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an Ermakov equation where real solutions must be chosen to make I Hermitian. oo 
is frequently rescaled to unity by a scale transformation of 6 [8J. Other common and 
convenient choice, which we shall adopt here, is uq = u(0). The eigenstates of I(t) 
become, with appropriate phase factors, solutions of the time-dependent Schrodinger 
equation, 

1/4 i 



V n (t,x) 



f mw \ 

V nk J (2 n n\b) 1 / 2 



exp 



x exp 



m / o 



X 



/ mu>o\ 1 / 2 x 
\~h~) b 



(38) 



(39) 



2h \b(t) b 2 

and form a complete basis to expand any time-dependent state, ip(x, t) = J2 n Cn^n(x, t), 
with the amplitudes c n constant. A method to achieve frictionless, population preserving 
processes is to leave u(t) undetermined first, and then set b so that 1(0) = Hq(0) and 
[I (t f) , H (t f)] = 0. This guarantees that the eigenstates of I and H are common at 
initial and finite times. We can do this by setting 

6(0) = 1,6(0) 



1,6 = 

b(t f ) = 7 = hM] 1/2 ,6(t/) = 0,b(tf) = 0, 



(40) 



and interpolating b(t) with some real function that satisfies these boundary condition. 
The simplest choice is a polynomial, 



b(t) 



(41) 



Once the aj are determined from (|40p . u(t) is obtained from the Ermakov equation ([3j 
and one gets directly a transitionless Hamiltonian Hu(t) = Ho(t) with a local, ordinary, 
harmonic potential, but note that uj(t) 2 may become negative for some time interval, 
making the potential an expulsive parabola [H[T8]. The II method is thus clearly distinct 
from from TT and implements a different Hamiltonian. Note also, by comparison of the 
coefficients, that the invariant operator / corresponding to Hjj is different from Htt, 
although they are both generalized harmonic oscillators. 



3.1. Physical realization 

The TT method only requires the time variation of a parabolic potential. Effective 
harmonic optical traps for neutral atoms may be formed by magnetic and/or optical 
means and their frequencies are routinely varied in time as part of many cold atom 
experiments. In magnetic traps, for example the frequency is modulated harmonically to 
look for collective excitation modes of a condensate [19] , and ramped down adiabatically 
to change its conditions (critical temperature, particle number, spatial extension) 
[SniElH], or as a preliminary step to superimpose an optical lattice [21 J . Some experiments 
involve both time-dependent magnetic and optical traps or antitraps [22J. Purely optical 
traps are also manipulated in time, e.g. for adiabatic cooling of single neutral atoms 
[23] . In particular laser beams detuned with respect to the atomic transition form 
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effective potentials for the ground state depending on Rabi frequency Vt and detuning 
A as Q 2 /4A by adiabatic elimination of the excited state, thus forming attractive or 
repulsive potentials. This effective interaction can be made time dependent by varying 
the laser intensity, the frequency, or both p], since the optical frequencies are many 
orders of magnitude larger than Rabi frequencies or detunings, and the changes will 
be slowly varying in the scale of optical periods. The intensity of a dipole trap can be 
changed by three or four orders of magnitude in 100 ns using acousto-optics or electro- 
optics modulators. To monitor the sign of the square frequencies, one can superimpose 
two dipole beams locked respectively on the blue and red side of the line. By controlling 
their relative intensity, one can shape the square frequencies and their signs at will. 

4. Discussion 

We have compared and distinguished two different methods: a "transitionless-tracking" 
(TT) algorithm, and an "inverse-invariant" (II) method, to achieve transitionless 
dynamics for a fast frequency change of a quantum harmonic oscillator. They imply 
different driving Hamiltonians. The one in the II method can be implemented for 
ultracold atoms or ions in the laboratory by varying the trap frequency in time along a 
certain trajectory, and a generalization to Bose Einstein condensates has been worked 
out [7j, but its extension to other potentials or systems may be difficult and remains 
an open question. By contrast, we have found some difficulties to realize the TT 
Hamiltonian for the harmonic oscillator, but the TT method has the advantage of being, 
at least formally, more generally applicable. The feasibility of the actual realization is 
quite another matter and has to be studied in each case. An example of application is 
provided in [3] for two-level systems. 
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Appendix A. Relation to the squeezing operator 

The evolution operator takes a particularly simple form when using the simplified case 
E n (t) = 0, so that H(t) = H\{t). Taking into account that [H\{t),H\{t')\ = we can 
write 



This may be evaluated explicitly with (fl~9l) fixing the time of the creation and 
annihilation operators to 0, 





(A.2) 
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which is a sqeezing operator with real argument r(t) = In ]. It is unitary with 

inverse [S^r)] -1 = S(—r). Using the relations 



a] + d t 



a] - d t 



uj(t) t „ 
<a + a ) 



\w(0) 



(A.3) 



and the formal properties of 5, see e.g. [21], it is easy to prove that 
S(r)a S(-r) = d t , 

S(r)dlS(-r) = at. (A.4) 

In fact any combination of powers of do and is mapped to the same combination of 
powers of d t and d\ by this unitary transformation. To show that |0 t ) = S^Oo) is indeed 
the vacuum at time t, note that 

a t |0 t ) = S{r)S{-r)aS{r)\0 o ) = S(r)d \0 ) = 0. (A.5) 

Similarly we note that, consistently, 

S(r)\n(0)) = S(r)^=(dl) n \0 ) = ^S(r)(dl) n S(-r)S(r)\0 ) 



a t 

n 



l)"|0 t ) = \n{t)). (A.6) 
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